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TECHNICAL NOTE 4346 

APPROXIMATE METHOD FOR CALCULATING MOTIONS IN ANGLES OF 
ATTACK AND SIDESLIP DUE TO STEP PITCHING- AND 
YAWING-MOMENT INPUTS DURING STEADY ROLL 
By Martin T. Maul and Teresa R. Brennan 

SUMMARY 


Hie simplified method of NACA Report 1344 for calculating motions 
in angles of attack and sideslip resulting from a trim angle of attack 
and steady rolling velocity has been extended in the present paper to 
the condition of pitching- and yawing -moment inputs. The resulting 
formulas are intended primarily for rolling-velocity conditions in 
which rolling divergence is not encountered. From the calculated angles 
of attack and sideslip, preliminary estimates of aircraft loads for 
design purposes can be made. Motions calculated with the simplified 
method are compared with exact solutions of the five-degree-of -freedom 
equations . 


INTRODUCTION 


Many high-performance aircraft encounter inertia coupling in 
rolling maneuvers. As a result this condition must be investigated 
during the design of new aircraft. A thorough analysis of rolling flight 
requires the solution of the nonlinear equations of motion, and analog 
or digi t al computers are usu ally utilized for this problem. However, 
for preliminary design purposes, simplified methods of analysis are 
desirable. 

A simplification used by many investigators in analyses of roll- 
coupled motions is that of linearizing the equations of motion by 
assuming rolling velocity to be constant. In early investigations of 
the stability of steady-rolling aircraft (refs. 1 and 2), the concept 
of critical ro lling velocities was defined. Steady-state solutions of 
these linearized equations also yield useful results as to the magni- 
tudes of the responses to be expected in rolling maneuvers and have 
been discussed in references 2 and 3. Transient solutions of the 
equations of motion for steady rolling can be obtained, but such 
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calculations are lengthy. Approximate expressions presented in refer- 
ence 4 permit simple calculation of transient responses of angles of 
attack and sideslip due to a trim angle of attack. 

In this paper the method of reference 4 is extended to the cal- 
culation of motions resulting from step pitching- and yawing-moment 
inputs. This method is applicable to inputs such as elevator deflec- 
tion, rudder deflection, aileron yaw, and yawing moment due to rolling 
velocity. Transient responses of angles -of attack and sideslip are 
presented for a range of rolling velocities for a current swept-wing 
fighter airplane assumed to be flying at a Mach number of 0.7 and an alti- 
tude of 32,000 feet. The method is intended primarily for conditions 
of rolling velocities up to the critical value. In the region of rolling 
divergence, nonlinearities in the equations of motion and airplane aero- 
dynamics are known to be significant, and the simplified, linearized 
results are not expected to be good approximations of the five-degree- 
of -freedom nonlinear results. 


SYMBOIS 


Airplane equations of motion are referenced to a system of prin- 
cipal body axes in this paper. 

A,B,C,D,E coefficients of equation (2) 

b wing span, ft 

c wing mean aerodynamic chord, ft 


C T lift coefficient, 

b qS 

C 7 rolling-moment coefficient, Roll ing moment 

1 qSb 

CL pitching-moment coefficient, Pitching agent 

m qSc 

C n yawing -moment coefficient, 

^ qSb 




side-force coefficient, 


Side force 
qS 


L X 


moment of inertia of airplane about X-axis, 


slug-ft 2 
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*Y 

h. 

Ka 


% 

Mz 

P 

Po 

q. 

r 

S 

t 

V 

a 

a Q 

/5a 


0 

^2 


o 

moment of inertia of airplane about Y-axis, slug-ft 

p 

moment of inertia of airplane about Z-axis, slug-ft 
amplitude coefficients (n = 0, 1, 2, and 3) 

airplane mass, slugs 
pitching moment, ft-lb 

yawing moment, ft-lb 

rolling Telocity, radians/ sec 

steady rolling velocity, radians/sec 

pitching velocity, radians/sec; dynamic pressure, lb/ sq ft 
yawing velocity, radians/sec 
wing area, sq ft 
time, sec 

airplane velocity, ft/sec 
angle of attack, radians 

initial angle of attack of airplane principal axis, radians 

incremental change in angle of attack, radians 

angle of sideslip, radians 

phase angle, radians 

real roots of characteristic equation 


a ± i£» complex roots of characteristic equation 

P 

“e = — ; 


2 
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•1 

ANALYSIS 


Equations of Motion 


If rolling velocity and forward speed are assumed constant and 
weight components are neglected, the airplane equations of motion refer- 
enced to principal axes are 


Pitching: 


Yawing: 


Side force: 


Normal force: 



( 1 ) 


The terms on the right-hand side of the equations are disturbances; 
the pitching-moment disturbance My results from elevator deflection, 

while the yawing-moment disturbance Mg may result from rudder deflection, 

N P 

aileron deflection (aileron yaw), or roll in g velocity — p Q . Airplane 

motions due to pitc hing - and yawing-moment inputs are the subject of this 
investigation. Motions due to the term UqPq were treated in refer- 
ence 4 and are omitted in this investigation. 

The characteristic equation for steady roll is given by the expan- 
sion of the dete rminant formed by the terms on the left-hand side of 
equations (l) and is 


AA^ + BA 5 + CA 2 + DA + E = 0 


( 2 ) 
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r 


where 


A = 1 


B _ i + ^ 

I z Ij aV mV 


C = - 


_ hz?l . h h 

mV Iy mV mV mV I2 



n r N r 2 Y B Y P M q, 
T + i + rr + p o + ”^ + -- 

% 1 Z 1 Y - L Z 


mV I z mV Iy 


Al r Mq\ 2_\^P + + 

^Z I yJ P ° % % h % P Z mV Iy I z mV I z mV 

!k i + bi ^ »r + ft i ^ + i tali') (; 

Iy mV mV mV I z mV mV Iy mV I z J\ 




E 





2 _ _ ^!£lP 

I Z P ° % I Z *Z mV I Y I Z mV 


Iy I z mV mV 


Y P p~X ~ ^Az ~ I x \_ 2 

mV mV\ I z /\ Iy ) 0 
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Approximate Expressions for Angles of Attack and Sideslip 

In reference 4, responses of angles of attack and sideslip due to 
the disturbance cypQ were obtained from equations (l) by the Laplace 

transform method ( ref . 5 ) • It was shown that the amplitude of the 
higher frequency mode appearing in the solutions was small and could be 
neglected. Thus, the solutions were simplified to a constant term plus 
either an oscillatory mode or two aperiodic modes, depending on the 
rolling velocity. In addition, the modal coefficients were simplified 
by omitting some terms of negligible value. The same procedure has 
been adopted In this investigation to derive approximate expressions 
for responses due to unit step inputs of My/ly and Mg^Ig. 


When the significant mode is a long-period oscillation, as for 
subsequent conditions In this paper with p Q = -1.0, -1-5, -2.5,> and 
-3.0 radians/sec, the response is written: 


P 

/ J Z 


P 

%/% 


V 1 * 

Ax 

%/% 


= + Kpe at cos(oat + 0) 


(3) 


When the significant modes are two exponentials, as for p Q = -1.86, 
-2.0, and -2.33 radians/sec in this paper, the response is written: 

p 


P 


A,t A p t 
= K Q + Kge L + K^e * 


(M 


V*J 
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After simplification, the amplitude coefficients are given by the 
following expressions. In these expressions, approximate values of the 
characteristic roots, as obtained by the method of the appendix of refer- 
ence 4, are used. 

Q 

In the expressions for — ~ — , 

w 

i + J5& + ( Iz " Ix V 2 

I Y mV I Y [ I Y fo 


K 


a 2\Ig 2 + H 2 

1 -l\Je/c 

where 


8 = tan' 


-1/ -HP - Ga \ 
\ Ha - (Sa>/ 


G = a) 2 - a 2 + ag - + KqE 


H = -2en> + cd 


V 


I Y mV , 


L = 4(3a 2 a> - ay 5 ) + 2Cu> 


K 2 " 


-h 2 + hi 

* mV j 

| + KqE 

A 1 (2CA 1 + D) 

1 


Kj - 


V + ^ ^ 


V 2 ( 2CA 2 


+ D1 
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In the expressions for 


M. 


■y/V 


% = 


Px - b) ^ . *r 
l Iz / mV I z. 


E 


K = 

d / e/c 


9 = tan’ 


-l/ -So - Ga \ 
l Ha - Qcd/ 


where 


G = p Q a 


H = p 0 cn 


- 

1 - 1 

(h - Vi 

-L 

Iz /J 

1 - 

r h - h\ 




V 


L = 4(3a^o - 0,3) + 20» 


Kg = 


Vi 

'*5^' 

1 

H 

1 i 

+ K q E 

V 2 ®! + D ) 



Pq A 2 

1 - 


+ KqE 

a 2 ( 

2CA 2 + 
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In the expressions for ; — , 

W 
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In. the expressions for 


Ax 


%/V 




N r Yp N p /I X - 2 

l Irr rO 


I z mV Ig 


E 


K - 2|G * 

K 1 


Lp/C 


6 = tan' 


•1 /-HCP - Ga \ 
\ Ha - Gen/ 


where 


G = 


2 2 
a - CD - a 


: [3 N. 


+ ^ + 


mV I z 


PP V 

H = 2an> - cd + -£ 
\mV In , 


% = 


L = 4(3a^D - or*) + 200 

^i 2 - A l(^ + ^) + ¥ X> E 
A-l^CA! + D) 


‘P t V 


K 3 


*2 ’ M^r + T^J + * 0 ® 

A^ ^2CAg + D ^ 
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DISCUSSION 


Airplane Stability in Steady Roll 

The stability of an airplane in steady roll is determined from the 
stability equation, equation (2), and has been treated generally in ref- 
erence 1 and, for the airplane of this investigation, in reference 4. 

The airplane physical and aerodynamic characteristics and flight condi- 
tions used in this paper are the same as in reference 4 and are given in 
table I. Exact and approximate roots of the characteristic equation 
were published in reference 4 and are repeated in table II. The rolling 
divergence boundaries, as developed in reference 1, are presented in 
figure 1 for this airplane. For this flight condition there is no diver- 
gent mode for any rolling velocity since the divergence boundary iB not 
intersected. At low roll rates (-1.0 and -1.5 radians/sec) the airplane 
has two stable osci lla tory modes. At roll rates for which the airplane 
is near the divergence boundary (-1.86, -2.0, and -2.33 radians/Bec), 
there is a stable osci llat ory mode and two convergent aperiodic modes, 
one of which is approaching neutral stability. At the higher roll rates, 
-2.5 and -3.0 radians/sec, the airplane again has two stable oscillatory 
modes. 


Sample Responses From Simplified Method 

Angle-of -attach and sideslip responses to unit step inputs of My/Iy 
and M z /I z have been calculated from the simplified expressions for 

a range of rolling velocities from -1.0 to -3-0 radians/sec and are shown 
in figures 2 and 3 along with the exact steady-rolling responses. These 
motions are for the airplane and flight condition given in table I. The 
good agreement between the simplified and exact motions justifies the 
use of the simplified expressions for steady-rolling calculations. 

The steady-state values of the responses shown in figures 2 and 3 
are given in table III. Since rolling maneuvers are generally of short 
duration, the steady- state condition is seldom attained. 

Note that for the larger roll rates of figures 2 and 3 only a 
portion of the responses shown might be of practical interest. For 
example, at a roll rate of -3.0 radians/sec a 3^0° roll would be performed 
in slightly over 2 seconds, whereas the responses are shown for a dura- 
tion of 8 seconds. 


Comparison With Five-Degree-of -Freedom Responses 

The nonlinear five-degree-of -freedom equations of motion are used 
universally for predicting airplane motions in large-anplitude rolling 
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maneuvers. In order to determine whether the simplified results are 
good approximations of the airplane motions experienced in practical 
rolling maneuvers, a comparison was made of the simplified and five- 
degree-of-freedom solutions. The airplane equations of motion with con- 
stant forward speed were programed on an analog computer and responses 
were obtained for a step aileron deflection in combination with (l) a 
step rudder deflection, (2) a step elevator deflection, and (3) combi n ed 
step elevator and rudder deflections. At the time of these calculations, 
more complete aerodynamic data which were available for this airplane 
(table IV } were used. The flight conditions and inertias were the 
same as for the calculations for figures 2 and 3* and the angle of attack 
of the principal axis was assumed zero. The results of these calcula- 
tions are presented in figure 4 with corresponding motions from the 
simplified method. 

The motions in angles of attack and sideslip are presented in fig- 
ure 4(a) for a step yawing -moment input - 0 . 2092 , equivalent 

to 5° of rudder deflection. The variation of ro lli ng velocity with 
time from the nonlinear equations is presented, from which trial values 
are selected for use with the simplified method. The values chosen were 
-2-5 radians/sec, the average value during a 360 ° roll, -3*0 radians/sec, 
and -3*2 radians/sec, the steady-state value. During the latter part 
of the motions the simplified method with a roll rate of -3-0 radians/sec 
provided the best agreement with the nonlinear motions. For the first 
second of the response, better agreement in the 0 motions is obtained 
when a roll rate of -1.5 radians/sec is used. This is the average value 
during the first second of the nonlinear rolling-velocity response. 

There is poor agreement initially in the angle-of -attack response because 
the longitudinal mode Is neglected in the simplified method. 

The response to a step pitching -moment input My/ly of 0.5879* 

equivalent to about -3° of elevator deflection, is shown in figure 4(b) 
for the nonlinear system and for the simplified method for roll rates 
of -1.5* -2.5* -2-7* and -3-0 radians/sec. The roll-rate response of the 
nonlinear system is not shown but is similar to the curve in figure 4(a) 
except that the steady- state value Is -2.7 radians/sec. For the first 
second the simplified results for a roll rate of -1.5 radians/sec, the 
average value, provide the best agreement with the nonlinear response. 

For the remainder of the time history, simplified results for the steady- 
state rolling velocity, -2-7 radians/sec, are good approximations of 
the nonlinear response. 

The results for combined pitching- and yawing-moment inputs are 
shown in figure 4(c). The steady-state rolling velocity is 
-2.6 radians/sec, and the nonlinear response is approximated in two 
steps as were the motions due to the individual inputs. 
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Suggested Procedure for Use of Approximate Method 

In airplane-design calculations, the simplified methods of this 
paper may "be employed as follows: Determine the max imum roll rate of 

interest and roll-rate response time for a specific flight condition 
from one-degree-of -freedom roll calculations and experience with s imilar 
aircraft. When the maximum roll rate is established, select intermediate 
values to cover the rolling-velocity range of the airplane. 

For a selected roll rate, plot the one-degree-of-freedom response 
( see sketch) : 


A 



For times up to t-^, compute Ax and (3 responses by using the 

average value of rolling velocity during this time. Compute the Ax 
and p responses for the steady-state roll rate as well for times 
beyond t-^. From these responses, maximum values of Ax and 3 for 
the maneuver are determined, with which preliminary load calculations 
may be made. A tabulation or plot of these data for the range of roll 
rates will disclose any serious conditions. Needless to say, critical 
conditions should be examined further by five-degree-of- freedom studies. 

A word of caution is advanced regarding the use of one-degree-of- 
freedom roll calculations. In roll-coupled maneuvers, the one-degree- 
of-freedom roll equation will not generally provide an accurate rela- 
tionship between aileron deflection and rolling velocity because of the 
effect of sideslip in damping or augmenting the roll. The use of the 
one-degree-of-freedom equation is intended only for obtaining estimates 
of response time and selecting a range of roll rates. Steady-state 
solutions of the five-degree-of-freedom equations are required to relate 
accurately aileron deflection and rolling velocity. 
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Angle of Attack of Principal Axis 

Angle of attack coupled with roll rate. is a major input to the 
equations for steady roll but was neglected in this paper in order to 
isolate the effects of pitching- and yawing-moment- inputs. The effect 
of angle of attack was treated in reference ^ and the results can be' 
algebraically added to the results of this paper for calculations in 
which both types of inputs are considered. 


Recovery Motions 

Recovery is the portion of the rolling maneuver in which the control 
is reversed to stop the roll., and the motions are defined by the non- 
linear equations of motion. In many cases the maximum excursions in 
angles of attack and sideslip (and, consequently, loads) occur during 
this t ime . For approximating recovery motions the classical three- 
degree-of-freedom lateral equations and two-degree-of -freedom longitu- 
dinal equations with initial conditions and approximations for the 
coupling terms were used. These equations are given in the appendix. 

Sanq)le recovery motions from the five-degree-of -freedom equations 
and two sets of approximate motions are presented in figure 5 * These 
recovery motions axe the latter part of a 360° roll after the aileron 
has been neutralized to stop the roll. Motions in angles of attack and 
sideslip, as calculated from the linear, uncoupled, two-degree -of -freedom 
longitudinal and three-degree-of -freedom lateral equations, are shown 
and differ in magnitude and phase with the nonlinear responses. When 
approximations to the coupling terms are included, the agreement with 
the nonlinear response is better. The remaining phase difference is 
attributed primarily to the 50 deg/sec rate at which the aileron was 
taken off In the analog calculations. When these approximate equations 
axe used for recovery calculations, occasional checks with the five- 
degree -of -freedom equations should be made. 

A dangerous recovery situation, for which the linear equations are 
invalid, is possible In flight. After the ailerons are neutralized for 
recovery from a rapid roll, it is possible for rolling velocity to be 
maintained at a large value instead of being reduced to zero. The air- 
plane enters an autorotational state, a steady-state condition of large 
angular end linear velocities. Five-degree-of-freedom solutions are 
required to predict this condition. 


Aircraft Loads 

In many rolling maneuvers the critical factor is the load encoun- 
tered. In pure rolls the vertical-tail load is usually the one of 
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concern. In a roll-coupled maneuver in which large excursions in angles 
of attach and sideslip occur, structural failures of the wing, horizon- 
tal tail, or vertical tail are possible. Preliminary estimates of sur- 
face loads can be made by using the calculated angles of attack and 
sideslip. 


CONCLUDING REMARKS 


The simplified method of MCA Report 1344 for calculating motions 
in angles of attack and sideslip in roll-coupled maneuvers has been 
extended in the present paper to the condition of pitching- and yawing-" 
moment inputs. The resulting formulas are intended primarily for rolling - 
velocity conditions in which rolling divergence is not encountered. 

Prom the calculated angles of attack and sideslip, preliminary estimates 
of aircraft loads for design purposes can be made. Motions calculated 
with the simplified method are compared with exact solutions of the five- 
degree-of -freedom equations . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 5, 1958 . 
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APPENDIX 

ROLL-RECOVERY EQUATIONS 


The five-degree-of -freedom, equations were separated into the three- 
degree-of -freedom, lateral and two-degree -of -freedom longitudinal equa- 
tions for approximating roll-recovery motions. As the coup lin g terms 
were known to have significant effects on recovery motions, they were 
included in simplified form. 

In the lateral equations the coupling terms appear as variables. 
These terms were linearized by assuming pitching velocity and angle of 
attack to be constant. After the Laplace transform is applied, the 
lateral equations are written as 


l p\_ fh: “ • L z'\_ _ "P 


r ■ w p ' i i 


x 


lo r - — P = P D 

ix 


fJ x- 




-a^ + r + (A - Jjp = P Q 


where is the yawing-moment input applied in this investigation 

and the subscript o denotes values of the variables at the start of 
recovery . 

In the longitu dina l equations the coupling terms appear as inputs 
since they consist of lateral variables only. The inputs represented 
by these terms are large at zero time and decrease to zero as the rolling 
velocity decays. Hence, these terms were approximated by assuming f3 
anri r to be constant and p to decay as given by the one-degree-of- 
freedom roll response. 

After the Laplace transform is applied, the longitudinal equations 
are written as 
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/, “A Hx 


H 

1 

Po r o 


A + ^ + l 

V ** / 

(A - A r ) 





-q + 

\ A + 5vy a " 0 

( x ■ 



where M^/l^ is the pitching -moment input applied in this investigation 
and Aj. is the one-degree -of- freedom characteristic rolling root. 
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TABLE I 

MASS AND AERODYNAMIC CHARACTERISTICS OF AIRPLANE 


m, slugs . . . 

l x , slug-ft 2 . 

2 

l y, slug-ft 

l z , slug-ft 2 . 

Ct , per radian 
a 



per radian 


Cjjj , per radian 
0 . 


C v , per radian 

C^, per radian 

C^ , per radian 
LL r 

V, ft/sec . . . 
S ; aq ft . . . . 


c, ft . • 
t, ft . . 
q, lb/sq ft 


745 

10,976 

57,100 

64,975 

5-85 

-0.36 

-5-5 

-0.28 

0.057 

-0.095 

691 

577 

11.3 

36.6 


197 





r 


TABLE II 

ROOTS OF THE CHARACTERISTIC EQUATION 


Exact 

Approximate 

- 0.488 ± 2.501 

-0.0729 ± 1.541 

- . 


-0.562 ± 2.891 

-O.199 ± 0.9421 

-0.552 ± 5-101 

-0.209 ± 0.8761 

- 0.537 ± 3-331 

- 0.224 ± 0.4851 

- 0.540 ± 5.421 

-0.221 ± 0.4691 

-O.527 ± 5.661 

-0.522, - 0.145 

- 0.532 ± 3-711 

-0.511, - 0.146 

- 0.324 ± 5-791 

-0.455, -0.020 

-0.529 ± 3-851 

-0.444, -0.020 

- 0.518 ± 4 .o 8 i 

-0.574, -o.m 

-0.522 ± 4.151 

-0.565, -0.113 

-0.516 ± 4.241 

-0.245 ± 0.2581 

-0.520 ± 4.291 

- 0.241 ± 0.2551 

-0.511 ± 4-701 

-0.250 ± 0.7601 

-0.512 ± 4.801 

- 0.248 ± 0.7431 


9^ HI TOVH 
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TABLE III 


STEADI-STATE VALUES OF MOTIONS PRESENTED IN FIGURES 2 AND 3 


radians 
V sec 

P __ _2 

2 

P 2 

Aa 2 

“z/V "" 


LL J v 

V 1 ! 

•yv 

-1 

-0.58 

-0.058 

-0.063 

0.213 

-1.5 

-1.07 

-.216 

-.236 

.249 

-1.86 

-3.57 

-1.35 

-1.47 

-.106 

-2.0 

.-13-34 

-7.03 

-7.65 

-3.45 

-2-33 

-.552 

-1.52 

-1.65 

-2.08 

-2-5 

.171 

-.503 

-.547 

-.898 

-3-0 

.21 

-.097 

-.106 

OJ 

00 

CJ 

1 
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FIVE-DEGREE-OF -FREEDOM PROBLEMS 

4.15 

-O.4o6 

-4.5 

-0.63 

0.0515 

- 0.175 

-O.325 

-0.03 

0.5879 


-0.2092 




V 

radians/sec 
O -1.0 
V -1.5 
□ - 1.86 
A -2.0 
O -2.33 
k -2.5 
zd -3.0 




(a-) p Q = -1.0 radian/sec. 


Time histories of ■ — and — for exact and simplified 

solutions of the four-degree-of -freedom equations. 
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0123 45678 

t, sec 


(b) P Q = -1*5 radians/sec. 


Figure 2.- Continued. 
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012345678 


t, sec 

( c ) p Q = -1.86 radians/sec. 


Figure 2.- Continued. 
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t, sec 


( d > P 0 


-2.0 radians/sec. 


Figure 2.- Continued. 
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■fc,- sec 


(e) p Q = -2-33 radians/sec. 


Figure 2.- Continued. 
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0 1 2 3 4 5 6 7 8 

t, seo 

(f) p Q = -2.5 radians/sec. 


Figure 2.- Continued. 
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t, BSC 

(g) P 0 = -3-0 radians/sec. 


Figure 2.- Concluded. 
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t, sec 

(c) p Q = -1.86 radians/sec. 


Figure 3-- Continued. 
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t, sec 


U) P 0 


-2.0 radlans/sec. 


Figure 5-- Continued. 
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0 l 2 3 k __ 5 6 7 8 

t, sec 

(e) P 0 = “2.33 radians/sec. 


Figure 3*~ Continued. 


NACA TN h3 46 


37 



t, sec 

(f) P D = -2.5 ra&ians/sec. 


Figure 3 .- Continued. 
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Figure 3*- Concluded. 
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(c) Combined yawing- and pitching-moment disturbances, — = O.5879, 
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— = = -0.2092. 
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Figure h.- Concluded. 
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Figure 5 .- Motions in angles of attack and sideslip during roll recovery. 
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